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Abstract 

In the low-density surface region of superfluid 4 He, the atoms are far apart 
and collisions can be ignored. The only effect of the interactions is from the 
long-range attractive Hartree potential produced by the distant high-density 
bulk liquid. As a result, at T = 0, all the atoms occupy the same single- 
particle state in the low-density tail. Striking numerical evidence for this 100% 
surface BEC was given by Pandharipande and coworkers in 1988. We derive a 
generalized Gross-Pitaevskii equation for the inhomogeneous condensate wave 
function 3>(z) in the low-density region valid at all temperatures. The overall 
amplitude of 3>(z) is fixed by the bulk liquid, which ensures that it vanishes 
everywhere at the bulk transition temperature. 
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I. INTRODUCTION 



Our starting point PJ is very simple, namely that at the free surface of superffuid 4 He, the 
local density of atoms n(r) becomes small. Since the system is like a dilute inhomogeneous 
Bose gas of 4 He atoms, the effect of inter-atomic collisions can be ignored. The dominant 
effect of the interactions is from the long-range attractive Hartree potential produced by the 
bulk liquid. The final result is that when n(r) is less than 10% of the bulk density at 
T = one has almost complete Bose-Einstein condensation (see Fig. 1 of Ref. 0). That 
is to say, all the He atoms in the low-density tail are in the identical single-particle state. 
This unusual surface Bose condensate should exhibit striking features when probed in an 
appropriate way (this "surface" condensate is not a 2D phenomenon, since the low-density 
region is spread over a large distance from the free surface). 

A key difference from the recently observed trapped Bose-condensed atomic gases is 
that through the long-range attractive tail of the interatomic potential, the bulk liquid 
determines the surface density profile and the amplitude of the Bose order parameter. The 
surface condensate is "driven" by the bulk liquid, which ensures that it vanishes in both 
the surface and bulk regions at the same bulk transition temperature. In this sense, the 
surface Bose condensate might be viewed as a realization of the "atom laser" which is being 
currently discussed in connection with trapped atomic gases. 

II. CONDENSATE WAVEFUNCTION IN SURFACE REGION 

Separating out the condensate contribution in the quantum field operators in the usual 
way, we have 

tp(r) = $(r) + ip(r) 

^t(r) = $*(r)+^(r), (1) 

where the condensate wavefunction $(r) = (^(r)) is the order parameter of a Bose- 
condensed system. The total density is given by 
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(^'r)^(r)) = n(r) = n c (r) + n(r), (2) 

where n c (r) = |$(r)| 2 is the local condensate and h(r) = (t{)<r)ip(r)) is the local density 
of 4 He atoms out of the condensate. The equation of motion for ip(r) leads to an exact 
equation for $(r) given by [[|.[IJ 

$(r) + J dr'v(r - r')(n(r')^(r)) = , (3) 

where v(r — r') is the bare He-He interaction. 

In a uniform Bose fluid, (^|) determines the bulk chemical potential /x in terms of $. In 
a system with a free surface in thermal equilibrium, one can find the asymptotic spatial 
dependence of $(r) where the density n(r) is small, since then (Q) reduces to 

$(r) = 0, (4) 

where 

Uij(r) = y dr'v(r — r')n(r') (5) 

is the Hartree potential at a point in the low-density region due to contributions from r' 
in the bulk region. Since |r — r'| is large, only the long-range attractive part of the He-He 
potential is relevant in (||). Contributions when r' is close to r are negligible since the density 
is so small (for further discussion, see Q). 

Far enough away from a planer surface centered at z — 0, Vjj(r) in (||) does not depend on 
the detailed form of the surface profile n(r). To leading order, we have Vh(z) = —nsa/z 3 , 
with a being determined by v(r). This allows us to solve (Q) for the asymptotic spatial 
dependence 

$( 2 ) = Ce -( Az+B / z2 \ (6) 

where A = \J 2n ^^ — lA _1 and B = mpBa/2h 2 A. While we claim this gives the exact 
asymptotic behaviour, we cannot calculate the amplitude C in (^) since this is determined 
by the bulk liquid. 
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The preceeding analysis leading to @ is valid at all temperatures below the bulk su- 
perfluid transition. We note that (Q) relates $(r) in the low-density region to the density 
profile n(r) in the bulk region, for which we do not have a microscopic theory. However, 
T = variational Monte Carlo calculations |2j carried out some years ago show that in 
the low-density region, n(r) approaches n c (r) = |$(r)| 2 and becomes identical with it when 
n(r) < dins. Thus, in this region (see Fig. 1 of Ref. Q), we have complete 100% BEC, 
with the local density given by the square of (g). 

While (U) is still valid, the amplitude C will decrease towards zero as we approach the 
bulk superfluid transition due to depletion of the condensate. At finite temperatures, the 
non-condensate local density n(r) will asymptotically approach the finite (but very small) 
uniform vapour density appropriate to a gas-liquid system in thermal equilibrium (This 
vapour is always a classical dilute gas since 3> ksT). 

As one gets near the interface, of course, n(r) will increasingly deviate from n c (r). We 
cannot solve (|3]) in this region. However, we call attention to the fact that the absolute 
value of the condensate density is expected to peak at the interface According to the 

numerical results of Ref. [JJ, one finds n c (r) ~ 0.2ng when n(r) ~ 0.5ns- This peak in the 
condensate density at the interface should be looked for experimentally. 



III. CONCLUSIONS 

We have argued that at T = 0, all the atoms are in the same single-particle quantum 
state in the low-density tail of superfluid 4 He. This inhomogeneous Bose-condensed gas is 
"controlled" by the bulk liquid through the long-range attractive tail of the He- He interaction 
potential. 

A time-dependent perturbation v ex (r,t) of the bulk liquid will drive this surface conden- 
sate coherently. Using the same physical arguments given in deriving (^), the equation of 
motion for $(r, t) = (tp(r)) t is 



d${r,t) 
dt 
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+ v ex (r,t) +v H (r,t) 



$(r,t) (7) 



with the time-dependent Hartree potential given by 

v H (r,t) = fdr'v(r-r'){n(r%. (8) 

As with (U), (|7|) is only valid for r in the low-density region and (Q) only includes the 
contributions from r' in the high-density bulk region. Experiments to probe (and utilize) 
this surface Bose condensate are needed, especially in connection with the unique coherent 
features it should exhibit. 
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